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A  Landauer  formula  for  the  current  through  a  region  of  interacting  electrons  is  derived  using  the 
nonequilibrium  Keldysh  formalism.  The  case  of  proportionate  coupling  to  the  left  and  right  leads,  where 
the  formula  takes  an  especially  simple  form,  is  studied  in  more  detail.  Two  particular  examples  where 
interactions  give  rise  to  novel  effects  in  the  current  are  discussed:  In  the  Kondo  regime,  an  enhanced 
conductance  is  predicted,  while  a  suppressed  conductance  is  predicted  for  tunneling  through  a  quantum 
dot  in  the  fractional  quantum  Hall  regime. 


PACS  numbers:  72.IO.Bf.  72.l5.Qm.  73.aO.Gk.  73.50  Bk 

The  formulation  by  Landauer  III  and  Biittiker  [2]  of 
the  current  through  a  finite,  possibly  disordered  region  of 
noninteracting  electrons  has  tremendously  enhanced  the 
understanding  of  transport  in  mesoscopic  systems  (3l. 
The  Landauer  formula,  which  expresses  the  current  in 
terms  of  local  properties  of  the  finite  region  (such  as  the 
transmission  coefficient)  and  the  distribution  functions  in 
connected  reservoirs,  has  been  used  extensively  and  sue* 
cessfully  in  many  areas,  including  the  scaling  theory  of 
localization  [4,5],  universal  conductance  fluctuations  [6], 
Aharonov-Bohm  conductance  oscillations  [7],  the  integer 
quantum  Hall  effect  (8l  and  its  quenching  [9],  the  quant¬ 
ization  of  ballistic  conductance  [10],  and  recently  in  the 
field  of  quantum  dynamics  of  driven  systems  (quantum 
chaos)  [I  ll. 

While  both  the  derivation  of  the  Landauer  formula  for 
noninteracting  electrons  [3]  and  its  application  are  well 
established,  an  apt  formulation  of  the  current  when  in¬ 
teractions  between  electrons  are  involved  has  been  lack¬ 
ing.  In  view  of  the  recent  technological  progress  in 
confinement  of  electrons  into  small  regions,  where  the 
electron -electron  interactions  plays  a  major  role  in  the 
transport  [12],  it  is  quite  clear  that  a  Landauer-type  for¬ 
mula  for  the  transport  through  such  an  interacting  region 
is  highly  desirable.  Several  attempts  have  been  matte  to 
deal  with  special  cases  (13).  Langreth  [14]  was  able  to 
express  the  linear  conductance  through  a  single  site  with 
an  on-site  interaction  (the  Anderson  model)  at  zero  tem¬ 
perature  in  terms  of  phase  shifts  and  thus  relate  the  con¬ 
ductance  to  a  scattering  matrix.  Apel  and  Rice  [15]  ap¬ 
proximated  the  interaction  in  one  dimension  by  the  values 
of  the  momentum  transfer  Sq  at  Sq  "0  and  Sq  ~2k/  and 
were  able  to  derive  a  Landauer-like  formula.  Unfor¬ 
tunately,  this  approximation  is  unsuitable  for  electrons 
confined  into  a  small  region.  More  recently,  Henhfield, 
Davis,  and  Wilkins  [16]  have  been  able  to  derive  a  for¬ 
mula  for  the  current  in  the  Anderson  model. 

In  this  Letter  (a)  we  derive  an  exact  formula  for  the 
current  through  a  region  of  interacting  cLctrons  coupled 
to  two  multichannel  leads  where  the  electrons  are  not  in¬ 


teracting.  The  formula  we  derive  [Eq.  (6)]  expresses  the 
current,  as  in  the  noninteracting  case,  in  terms  of  the  Fer¬ 
mi  functions  in  the  leads  and  local  properties  of  the  in¬ 
teracting  region,  (b)  We  show  how  the  noninteracting 
case  and  the  results  of  Langreth  and  of  Henhfield,  Davis, 
and  Wilkins  follow  as  special  cases,  (c)  The  current  will 
be  written  in  a  particularly  simple  form  for  the  case  of  a 
constant  asymmetry  factor  relating  the  coupling  to  the 
left  lead  to  the  coupling  to  the  right  lead.  This  case  wifi 
be  presented  and  investigated  in  some  detail,  (d)  Two  ex¬ 
amples  where  the  interactions  lead  to  interesting  results 
for  the  current  will  be  discussed:  transport  in  the  Kondo 
regime,  where  an  enhanced  conductance  has  been  predict¬ 
ed  [17],  and  tunneling  through  a  correlated  electronic 
state,  such  as  the  fractional  quantum  Hall  state,  where  a 
suppression  of  the  conductance  is  expected  [1 8]. 

Our  starting  point  is  the  Hamiltonian 

Z  wU*.+//«,({d.t];]dJ) 

*.«  €  l.k 

+  Z  (K*McU,+H.c.),  (I) 

m 

where  cfa  (c*«)  creates  (destroys)  an  electron  with 
momentum  k  in  channel  a  in  either  the  left  (£,)  or  the 
right  (R)  lead,  and  idil  and  id,}  form  a  complete,  ortho¬ 
normal  set  of  single-electron  creation  and  annihilation 
operators  in  the  interacting  region.  The  channel  index  in¬ 
cludes  spin  and  all  other  quantum  numbers  which,  in  ad¬ 
dition  to  k,  are  necessary  to  define  uniquely  a  state  in  the 
leads.  This  Hamiltonian  corresponds  to  an  experimental 
situation  (see,  e.g..  Refs.  (10)  and  [12])  where  two  metal¬ 
lic,  multichannel  leads  are  connected  to  the  mesoscopic 
system  under  study  (see  Fig.  I ).  Since  the  leads  are  me¬ 
tallic,  the  interaction  between  electrons  in  the  leads  and 
the  interaction  between  electrons  in  the  intermediate  re¬ 
gion  and  electrons  in  the  leads  are  strongly  screened  and 
can  be  neglected.  However,  in  the  presence  of  barriers 
between  the  leads  and  the  intermediate  regr",  electrons 
cannot  flow  freely  to  screen  the  interactions  in  the  inter¬ 
mediate  region.  Accordingly,  the  interactions  between 
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FIG.  I.  Schematic  diagram  of  the  experimental  configu¬ 
ration  for  which  an  interacting  Landauer  formula  for  the 
current  is  derived.  Two  leads,  characterized  by  chemical  poten¬ 
tials  tn  and  tm,  are  connected  to  a  mesoscopic  region  where 
electrons  may  interact.  If  pi.  >  tin.  an  electron  current  J  will 
flow  from  left  to  right. 

electrons  in  the  intermediate  region  need  to  be  treated 
dynamically  and  are  included  in  Eq.  (I).  This  treatment 
is  similar  to  the  one  leading  to  the  Anderson  model  of  a 
single  magnetic  impurity  in  a  metal  1 1 9], 

The  logic  of  our  approach  follows  that  used  in  Ref. 
[20]  for  the  one-dimensional  noninteracting  case.  The 
unperturbed  system  (taken  to  exist  at  (■-«)  consists 
of  three  uncoupled  regions:  a  left  lead  and  a  right  lead, 
both  described  by  the  first  term  in  (I),  and  an  interact¬ 
ing,  intermediate  region  described  by  the  second  term  in 
( I ).  Since  the  leads  are  not  coupled  at  /  ■■  —  «>,  each  one 
maintains  its  own  thermal  equilibrium  and  one  can  asso¬ 
ciate  chemical  potentials,  pi  and  pg,  with  the  left  and 
right  leads,  respectively.  As  the  coupling  turns  on  be¬ 
tween  the  intermediate  region  and  the  two  leads  [the  last 
term  in  (1)1.  then,  if  pL  >  Ptt>  an  electron  current  J  starts 
to  flow  from  the  left  lead  to  the  right  lead.  After  some 
time  the  system  achieves  a  steady  state.  Our  aim  is  to  re¬ 
late  the  steady  state  current  to  pi  and  fig,  or  equivalent¬ 
ly,  to  /t(e)  and  /*(<),  the  unperturbed  Fermi-Dirac  dis¬ 
tribution  functions  in  the  leads.  (We  assume  that  the 
reservoirs  are  large  enough  that  the  bulk  pi  and  fig  are 
not  perturbed  by  the  current  /.) 

To  this  end  we  write  the  current  between  the  inter¬ 
mediate  region  and  the  left  lead  as 

J-- 1 

rl  k,m€L 
m 

ft  i  “  2k 

19 

(2) 

The  first  line  in  (2)  can  be  easily  checked  by  writing  the 
continuity  equation  for  the  current  [20],  while  to  get  the 
second  line  we  have  used  the  definition  of  the  Keldysh 
Green  function  [20,21],  G,>«(/)|i/(c^«di,(/)).  In  the 
Keldysh  formalism,  since  the  Hamiltonian  describing  the 
leads  is  noninteracting,  one  has  the  Dyson  equations 


Gkm.n  (oj)  “2  Kla  (to  JG.f,  (a>  ) 

m 

_  (J) 

Gm<i,(at)  “J  V*m  „(w) 

m 

~gL.it(a*)G„<m(a/)l . 

where  G.<*(/)«i(d£d)t(f)>.  and  G„>„(i),  to  be  used  later, 
is  equal  to  — /(dnOJdi).  The  Green  functions  with  su¬ 
perscripts  t  and  i  are  the  time-ordered  and  the  anti-time- 
ordered  Green  functions,  respectively  [22],  and  the  Green 
functions  denoted  with  small  g  are  the  unperturbed  Green 
functions  (i.e..  in  the  uncoupled  system).  Using  the 
equalities  [22]  G  >(w)+G  <(a»)“G'(a>)  +  G'(w)  and 
G  >(a>)  —  G  <(t>)mmG'(a>)~G*((o),  where  G'  (G®)  are 
the  usual  retarded  (advanced)  Green  functions,  and  the 
relations 

ft£*«(®)“2*//i.(»)$(«-*i.) , 

(4) 

&«.*•(»)•  ~2*»ll  —  fL(a>)]8(at-  (kt)  , 
where  a  €  L,  we  find 

ft  «€t*' 

x(/-,.(e)[G;.«(f)-G;.«(f)]+G.<m(e)l , 

(5) 

where  pm(e)  is  the  density  of  states  in  channel  a  and 
V ..«(f )  equals  for  An  equivalent  formula 

can  be  derived  for  the  current  between  the  intermediate 
region  and  the  right  lead.  Since,  in  steady  state,  the 
current  is  uniform,  one  can  symmetrize  Eq.  ($),  and  us¬ 
ing  matrix  notation  for  the  level  indices  in  the  interacting 
region,  we  find 

J  -  /  dt<u{[fL{t)rL  -/*(<)  r*l(G'  -  G°)i 

+tr{(rt-r*)G<)).  (6) 

where  r£«,-2 zP.(f)F«J,(#)F^t(r),  with 
defined  similarly.  In  equilibrium,  /*,(*)  *■/*(*)  —/«q(<r), 
one  has  G  <  “  — /«,( G'- G*)  and  the  current  vanishes. 

Equation  (6)  is  the  central  result  of  this  work.  It 
expresses  the  current  through  the  interacting  region  in 
terms  of  the  distribution  functions  in  the  leads  and  local 
properties  of  the  intermediate  region,  such  as  the  occupa¬ 
tion  and  the  density  of  states.  (The  local  density  of  states 
is  proportional  to  the  diagonal  part  of  G'-G",  while  G  < 
is  a  product  of  the  density  of  states  and  the  occupation 
[23].)  Note  that  these  are  to  be  calculated  in  the  pres¬ 
ence  of  the  leads. 

For  the  noninteracting  case  one  can  write  down  Dyson 
equations  for  the  Green  functions  in  the  intermediate  re¬ 
gion  [20],  G  <  m,ift.(<)GTLQa+ifg(e)GT*Ga  and  Gr 
—  G'“  —  iG,(r1+r*)G",  which  enable  us  to  rewrite 
(6)  as 

J— *-  fdelf, .(()  -/*(«)!  tr)GT*GTM  (7) 
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Since  the  transmission  coefficient  from  left  to  right  is 
given  by  f«.«  52xL,  p'J2 p'J' with  a  6  R 
and  a  €  L,  Eq.  (7)  reduces  to 

J  -/#(f)]tr{tt*(f)l  .  (8) 

which  is  the  usual  two-terminal  Landauer  formula  for  the 
noninteracting  case  (24). 

The  conductance  formula.  Eq.  (6).  takes  an  especially 
simple  form  {25]  for  the  case  that  the  couplings  to  the 
leads  differ  only  by  a  constant  factor.  rL(e)  ■*xr*(r), 

/*tA<e)-/*(e)Jtr{r<G'-G')} 


at  zero  temperature.  At  finite  temperature,  or  at  finite 
voltage,  inelastic  processes  have  to  be  taken  into  account 
and  the  usual  Landauer  formula  (8)  breaks  down  for  an 
interacting  system. 

In  order  to  demonstrate  the  new  features  of  Eq.  (9).  we 
now  study  two  specific  examples.  The  first  is  transport 
through  a  quantum  dot  in  the  Coulomb  blockade  regime 
[1 21.  Recently,  the  Anderson  Hamiltonian  {19], 

H  ™  £eiac/9Cia  +  <oXd^d<,  +  (/n)iil 

k.a  9 

+  I  O'taCU.+  H.C.).  (||) 

k  €  I..K 


-/#(*)] ImltdrGl] ,  (9) 

where  r*I'tr'7(rt  +  r*).  While  this  innocent  looking 
formula  resembles  the  noninteracting  one.  it  should  be 
emphasized  that  even  though  there  is  a  single  integral 
over  energy  in  (9),  this  formula  includes,  by  means  of  the 
full  Green  function  G',  inelastic  processes,  spin  flips,  and 
even  processes  where  several  electrons  are  scattered. 

In  order  to  illustrate  how  the  additional  processes  due 
to  interactions  are  reflected  in  the  formula  for  the 
current,  we  use  the  usual  definition  of  the  self-energy  £', 
namely,  £'m(g') _l  —  (G') -  l.  with  I*  defined  similarly, 
to  write  G'-G*-GTG*.  where  £-£'-£'.  This  al¬ 
lows  us  to  rewrite  Eq.  (9)  as 

Jmjf  /l(«)  "  /«(<)Ur{rGTGl 
" -jj-  J* ddfL (e)  —/#(<)]  tr{GT*GTtI<f  'll .  (10) 

where  Eo.  the  self-energy  for  the  noninteracting  case,  is 
equal  to  —  /'(r^+r*).  Comparing  Eq.  (10)  to  the 
noninteracting  results,  Eqs.  (7)  and  (8).  we  see  that  in 
the  presence  of  interactions  the  current  cannot  in  general 
be  recast  in  terms  of  the  transmission  matrix,  and  that 
the  additional  processes  included  in  £  are  directly  respon¬ 
sible  for  the  deviation  of  the  exact  formula  for  the  in¬ 
teracting  case  from  the  usual  Landauer  formula  (Eq. 
(8)].  Note,  however,  that  at  zero  temperature  and  in 
linear  response  only  single-electron,  elastic  processes  are 
allowed  by  energy  conservation.  In  this  case  £(p) 
“£o(p )  and  Eq.  (10)  reduces  to  the  usual  Landauer  for¬ 
mula  (8)  (this  is  a  generalization  of  the  Langreth  (14]  re¬ 
sult  for  the  Anderson  model).  Thus  the  Landauer  formu¬ 
la  for  the  linear-response  conductance  holds  not  only  for 
the  noninteracting  case,  but  also  for  the  interacting  case 


has  been  employed  (26]  to  describe  transport  in  this  re¬ 
gime.  Equation  (II)  is  a  special  case  of  (I),  where  the 
different  channels  are  the  two  spin  directions,  and  the  in¬ 
teracting  region  is  a  single  site  with  an  on-site  Coulomb 
repulsion  V.  In  this  case  the  current  takes  the  form 


•f-f  1/ <MA(«) -/*(«)  lr„„(<) 


--lmG^(f) 

a 


(12) 

where  -(l/*)ImG£.(e)  is  just  the  local  density  of  states 
of  electrons  with  spin  a.  In  linear  response,  for  tempera¬ 
tures  larger  than  the  Kondo  temperature,  the  conduc¬ 
tance  will  consequently  exhibit  resonant  tunneling  peaks 
only  at  to  and  «o+(/.  which  correspond  to  resonances  in 
the  density  of  states  of  the  uncoupled  site  (26].  However, 
below  the  Kondo  temperature,  the  density  of  states  devel¬ 
ops  a  peak  at  the  Fermi  energy  (27],  for  *o  <  m  <  (q+U. 
Consequently.  Eq.  (12)  predicts  a  greatly  enhanced  con¬ 
ductance  over  this  entire  range,  in  agreement  with  earlier 
studies  (17].  Equation  (12).  which  has  been  indepen¬ 
dently  derived  by  Hershfield,  Davis,  and  Wilkins  (16]  for 
constant  T,  provides  a  framework  to  study  the  crossover 
from  high  to  low  temperatures,  and  to  calculate  the 
current  in  a  Kondo  system  out  of  equilibrium.  Detailed 
studies  will  be  presented  elsewhere  (28). 

A  second  example  where  nontrivial  effects  due  to  in¬ 
teractions  appear  in  the  conductance  is  the  case  of  tunnel¬ 
ing  through  a  quantum  dot  in  a  highly  correlated  state, 
such  as  a  fractional  quantum  Hall  state.  In  the  case 
where  the  coupling  to  the  leads  is  weak,  i.e.,  the  elastic 
broadening  of  the  levels  is  smaller  than  the  excitation  en¬ 
ergy,  the  correlated  eigenstates  are  only  weakly  perturbed 
by  the  leads.  When  the  temperature  is  larger  than  the 
elastic  broadening  one  can  rewrite  Eq.  (9)  for  the  linear- 
response  conductance  G  in  the  form 


G-- ^-2r,.m(£y 


■E^'ZiPi+Pj) 

i.J 


(£>-£,) 


(13) 


where  the  are  eigenstates,  with  energies  £,,  of  the  uncoupled  interacting  region,  and  Pi  is  the  equilibrium  probability 
of  state  fi.  Hence,  the  conductance  consists  of  thermally  broadened  resonant  tunneling  peaks  which  occur  whenever  the 
chemical  potential  n  suffices  to  add  another  electron  to  the  interacting  region.  For  noninteracting  electrons,  one  can 
choose  the  d’s  to  correspond  to  single-particle  eigenstates  of  the  uncoupled  system,  and  the  overlap  factor  in  each  term, 
^y|d,lr>><r.  ldm|r/).  is  trivially  0  or  I.  On  the  other  hand,  adding  an  electron  in  a  single-particle  state  to  a  correlated 
V-particle  eigenstate — which  cannot  be  written  as  a  Slater  determinant  of  single-particle  states — will  generally  not  pro- 
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duce  an  (A/  + 1  )-particlc  eigenstate,  and  consequently 
these  overlap  factors  can  reduce  the  conductance 
significantly.  In  particular,  the  conductance  through  a 
quantum  dot  in  the  fractional  quantum  Hall  regime  has 
been  studied  in  detail  in  Ref.  (18)  where  it  was  shown 
that  the  overlap  factors  reduce  the  conductance  peaks  by 
a  factor  of  I /N<*~t)n,  for  the  v—  i/p  Laughlin  state, 
where  /V  is  the  number  of  electrons.  Thus  Eq.  (13)  pre¬ 
dicts  a  strong  suppression  of  the  conductance  in  this  re¬ 
gime,  an  effect  of  considerable  experimental  significance. 
Moreover,  formula  (13)  provides  a  way  to  calculate  the 
conductance  in  more  complicated  fractional  quantum 
Hall  states  (such  as  v“  f  ),  where  even  more  interesting 
effects,  such  as  an  alternating  suppression  of  the  conduc¬ 
tance  peaks,  are  predicted  to  occur  (l 8]. 

To  conclude,  we  have  derived  a  Landauer-type  formula 
for  the  current  through  an  interacting  electron  region. 
This  formula  provides  a  new  framework  to  study  trans¬ 
port  in  mesoscopic  systems  where  interaction  effects  are 
important.  We  have  demonstrated  the  usefulness  of  this 
new  tool  for  interacting  systems  by  identifying  novel 
features  in  the  conductance  for  two  examples:  the  Kondo 
regime  and  the  fractional  quantum  Hall  regime.  Further 
investigations  are  currently  under  way. 
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We  discuss  the  effects  of  surfaces  on  bulk  plasma  oscillations  in  a  metallic  film  of 
thickness  L.  The  leading  corrections  to  the  plasma  frequencies  due  to  the  surfaces  are 
proportional  to  L~l.  The  frequencies  of  low-lying  long  wavelength  modes  are  given 
by  the  dispersion  relation  of  bulk  modes,  w  =  wp(l  +  a(q2/k2F) ],  where  the  allowed 
q- values  are  q^  —  (nr/T),  T  =  L  -  2 dt,,  and  d(,  is  a  complex  length  characteristic  of 
the  surface  (reminiscent  of  but  different  from  Feibelman’s  d±).  An  explicit  expression 
for  df,  is  derived.  Resonant  excitation  by  an  external  field,  is  calculated. 


1 


The  widespread  interest  in  nanostructures  (thin  films,  layer  structures,  quantum 
wires,  quantum  dots  etc.)  as  well  as  in  atomic  and  molecular  clusters,  in  which  at 
least  one  dimension  is  “mesos topic”  (10  A  -  1000  A),  has  led  us  to  re-examine  the 
effects  of  finite  size  on  charge  collective  modes  (plasmons).  In  an  infinite  system, 
these  modes  have  frequencies  given  by  the  bulk  plasmon  relation  [lj 


u(q)  =  wp  +  a 


(1) 


where  uv  is  the  classical  plasma  frequency,  =  iire2n/m,  with  n  the  electron  den¬ 
sity,  kp  the  Fermi  wavevector,  and  a  a  (complex)  number.  The  imaginary  part  of 
a  describes  the  damping  of  bulk  plasmons,  which  in  an  infinite  system  is  due  to  ex¬ 
citations  of  multiple  electron-hole  pairs.  In  a  finite  system,  we  expect  the  following 
qualitative  changes  in  the  plasmon  spectrum  due  to  the  presence  of  surfaces.  First, 
only  a  discrete  set  of  frequencies  will  be  allowed.  Second,  there  will  be  additional 
damping  due  to  the  decay  of  plasmons  into  (single)  electron-hole  pairs.  In  this  Let¬ 
ter,  we  will  show  how  both  the  real  and  imaginary  parts  of  the  allowed  frequencies 
are  completely  determined  by  a  single  (complex)  length.  Specifically,  we  consider  a 
thin  slab  with  a  positive  jellium  charge  background,  n+  =  ft  between  2  =  0  and  z  =  L 
and  infinite  in  the  x  and  y  directions,  and  a  neutralizing  electron  liquid  described  by 
the  equilibrium  density  distribution  n_(z)  (See  Fig.  1).  We  shall  discuss  bulk-like 
collective  modes  of  the  form 


n(r,()=n(2)e-^*”,  (2) 

where  the  eigenfrequencies  u/n)  are  near  the  classical  bulk  plasma  frequency  associated 
with  the  density  n,  o/p  =  (4 Tre2n/m)1/2. 

This  problem  has  been  previously  addressed  experimentally  [2,3]  and  theoretically 
[4]  for  modes  of  the  form  (2)  and  of  the  more  general  form 

n(r,t)  =  n(z)e,<p*-'*>.  (3) 


In  the  present  Letter,  limited  to  the  case  p  =  0,  we  derive  the  following  new  result: 
the  (complex)  eigenfrequencies  u/n),  near  u/p,  are  determined  by  the  bulk  plasmon 
dispersion  Eq.  (1)  for  small  q,  and  by  an  effective,  complex,  slab  thickness 

T  -  L  —  2dj,  (4) 

where  db,  precisely  defined  below,  is  an  effective,  complex,  surface  thickness.  The 
thickness  db  is  reminiscent  of  the  parameter  d±  introduced  by  Lang  and  Kohn  [5]  for 
u>  =  0  and  generalized  for  arbitrary  u>  by  Feibelman  [6]  but  is  not  the  same  quantity. 
The  allowed  (complex)  interior  wavenumbers,  q^n\  are  given  by 

q^T  =  nir  (n  =  1,2,...)  (5) 


and  the  corresponding  (complex)  eigenfrequencies  are  given  by  Eq.  (1): 


2 


(6) 


J"  =^(9ln>). 

We  now  derive  the  above-stated  results.  We  take  the  equilibrium  electron  density, 
n_(z),  to  vanish  for  z  <  -c  and  z  >  L  +  c,  and  introduce  an  auxiliary  thickness 
[7]  a,  large  compared  to  the  surface  thickness  but  much  smaller  than  L.  We  then 
distinguish  two  surface  regions,  — c  <  z  <  a  and  L  -  a  <  z  <  L+c,  and  the  bulk 
region  a  <  z  <  L  —  a,  in  which  n_(z)  =  n  (see  Fig.  1). 

We  consider  a  collective  mode  of  frequency  u >,  with  self-consistent  electric  field, 
£(z),  parallel  to  z.  This  field  induces  a  current  density, 

j(z)  =  £‘dz'c(z,z')E[z'),  (7) 

where  <7(z,  z')  =  f  dx1  dy1  <r,,(x  —  x\y  —  y’\  z,  z');  att  is  the  zz-component  of  the  non¬ 
local  conductivity  tensor.  (For  ease  of  notation,  the  time  dependence,  e~‘wt,  and  the 
functional  dependence  of  fields,  currents  etc.  on  u  will  generally  not  be  explicitly 
written  out.)  The  quantity  <r(z,z')  is  short-ranged  in  the  difference  variable  (z  —  z')\ 
also,  when  both  z  and  z'  are  in  the  bulk  region, 

<t(z, z')  =  a{z  -  z'),  a<  z,z‘  <  L-a  (8) 

equal  to  the  bulk  conductivity  corresponding  to  n  and  u>.  Of  course  <r(z,  z')  vanishes 
when  either  z  or  z/  are  outside  the  slab  (— c,  L  +  c),  as  do  E(z)  (due  to  charge 
neutrality)  and  j(z). 

Denoting  the  induced  charge  density  by  n(z),  the  continuity  equation  and  Gauss’ 
law  are 

iuen(z)  +  =  0  (9) 

dz 

and 

d-y?-  =  _4jren(z)  (10) 

dz 

(where  e  has  been  taken  as  positive).  Since  at  z  =  -c  both  j(z)  and  E(z)  vanish,  (9) 
and  (10)  imply  that,  for  all  z, 

j(z)  =  ~E(z).  (11) 

We  can  eliminate  j(z)  from  (7)  and  (11)  which  results  in  the  integral  equation 

jL_*‘ dz' a(z,z')EW).  (12) 

We  are  looking  for  solutions,  which,  in  the  bulk,  behave  as  e±,qx,  with  q  <C  kp. 
In  this  limit  of  long  wavelengths  in  the  bulk,  we  can,  with  Eq.  (8)  in  Eq.  (12)  obtain 
the  frequencies  of  these  solutions  by  expanding  E(z')  about  z  when  z  and  z'  axe  in 
the  bulk.  The  result  is 


3 


(13) 


where 


ILJ 

4tt 


=  <ToM  -  v^')92> 


t  itjJ* 

<Jo(u)  =  J  dz  (j{z  -  Z\u)  =  J-2-, 
<r2( u/)  =  J dz’(z'  —  z)2<r{z  -  z';u/). 


(14) 

(15) 


Equation  (13)  is  the  bulk  plasmon  dispersion  relation  of  an  infinite  svstem  of  conduc¬ 
tivity  <r(z  —  z').  In  the  limit  L  >  a,  the  frequencies  of  the  lowest- lying  eigenmodes 
of  Eq.  (12)  are  very  near  uv  and  the  shape  of  E(z)  in  the  surface  region  tends  to  the 
limiting  shape 

E(z)-E(z)mE(r,  up)  (16) 


shown  schematically  in  Fig.  2.  We  arbitrarily  take  E(a)  =  1  and  write,  near  z  —  a 


E(2)  -  ~  (17) 

fl  — 

where  db  is  the,  as  yet  undetermined,  point  where  E(z)  extrapolates  to  0.  When  z  is 
in  the  bulk  the  solution  which  joins  (17)  near  z  =  a  is 

E(z)=  v-1  sinq(z-db).  (18) 

q(a  -  db ) 

Thus,  db  plays  the  role  of  a  (complex)  scattering  length. 

The  length  db  can  be  simply  expressed  in  terms  of  the  limiting  surface  density 
distribution,  n(z),  associated  with  collective  modes  with  q  — ►  0  and  u>  — ►  uy  From 
Eq.  (10)  and  the  form  of  E(z)  in  the  surface  region  we  see  that  n(z)  has  the  general 
form  shown  in  Fig.  2.  Thus,  for  z  near  a,  we  have 


— c 

,  dn(z') 


dz 1 


dz'2'  d*> 


J-c  dz  J-e 


00  dz’ 


dz' 


(19) 


The  upper  limits  oo  are  to  be  understood  as  >  k^x  but  <  L.  The  linear  form  of 
E{z)  near  a,  Eq.  (19),  extrapolates  to  sero  when  z  equals 


(20) 
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Thus,  db  is  the  center  of  mass  of  the  function  dn(z)/dz. 

The  eigenmodes  are  either  even  or  odd  about  L/2  so  that,  by  Eq.  (18),  the 
allowed  values  of  q  are 

qln)T  =  nx,  (n=  1,2,...)  (21) 

where  X  is  give  by  Eq.  (4).  The  corresponding  frequencies  are  given  by  Eq.  (1). 
Using  the  notation  db  =  db  +  idb  etc ,  we  have 


u> 


(n)" 


(22) 

(23) 


where  L  =  L  —  2dfb.  We  observe  that  the  real  parts  Eq.  (22)  of  the  frequencies 
correspond  to  bulk  plasmon  modes  with  the  boundary  conditions  E{d!b)  =  E(L—d'b)  = 
0.  The  imaginary  parts,  (23),  have  contributions  from  the  bulk  (oc  q")  and  from  the 
surface  (oc  (d"f X*)). 

We  have  solved  the  integral  equation  (12)  in  the  discrete  form 

iu  * 

—  Em  -  b  52  (24) 

n=0 

where  b  s  L/N  is  the  interval  length.  We  have  chosen  L  =  35.4  A,  N  =  248,  b  =  0.14 
A,  and  c  =  0.  For  the  sake  of  illustration,  we  have  chosen  the  following  values  for 
<rmTl.  In  the  “bulk”,  4  <  m  <  125,  we  have  chosen 

*Ql]  >  crm,m+3  =  O’Vn.m-a  — 

*Ql]  »  &m,m+ 3  =  &m,m— 3  =  ^ (25) 

and  symmetrically  for  m  >  125,  with  kp  —  1.75  x  10®  cm-1,  hvp  =  15  eV,  and 
a  =  0.57  +  t0.025  in  Eq.  (1)  (appropriate  [1]  for  A1  in  the  limit  of  q/kp  <  1).  The 
functions  n(u)  and  v  are  chosen  so  that  the  average  and  second  moment  of  <rmn  are 
<70(u;)  and  <72( up).  In  the  “surface”,  0  <  m  <  3,  we  have  chosen 


^m.m  —  t/i(w)  — 


IU)‘ 


4  iru>b  3  irkibP 


^m,m+ 1  —  ®m+I,m  —  — 


2wd 


2l7rAri.fr3 


[o2  - 
(o2  — 


<7'on  —  <7n(j  —  0 


m 


<7m,m  =  —  (*>(ta»)  +  6) 


m 


^mtm+l  —  —  .  {}V  4*  c) 

4 


m 


2  —  m+2,m  —  g  (*'*'  4"  c) 


m 


771+3  —  ^771  +  3,771  —  (i*7  4“  c) 


(26) 
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and  symmetrically  near  the  other  surfaces.  The  small  real  part  t  =  *>pa" /{fiTrkjrb3) 
was  added  to  the  conductivity  in  the  surface  regions  as  a  simple  model  of  damping 
in  the  surfaces  due  to  decay  into  electron-hole  pairs.  Apart  from  this  small  real 
part,  equation  (26)  is  a  simple  interpolation  between  aQn  =  <rnQ  =  0  and  the  bulk 
conductivity.  The  results  for  the  lowest  three  modes  are  plotted  in  Figs.  3  and  4,  and 
confirm  our  conclusions,  Eqs.  (5),  (6),  (22)  and  (23). 

How  do  our  results  relate  to  charge  density  waves  in  classical  Maxwell  theory? 
“Classical”  here  means  a  local  a (z,z') 

iuP 

■  <^ciw(2,  z')  =  (To 6(z  -  z)  =  -  *')»  (27) 

so  that  (12)  becomes 

|£W  =  ^£W.  (28) 

The  boundary  conditions  [8]  are,  from  charge  neutrality,  £(0)  =  E(L)  =  0.  Thus, 
classically  the  general  solution  is  u>  =  uv  and  E(z )  =  F(z ),  an  arbitrary  function  of  z 
satisfying  the  boundary  conditions  E( 0)  =  E(L)  =  0.  In  particular,  the  functions 

F„(*)  =  sin(g£i,z)  gJi  = —,  n  =  0,1,...  (29) 

£<">=u>p,  all  n,  (30) 


are  solutions.  Thus,  our  non-local  and  translationally  non-invariant  <r(z,  z')  leads  to: 
(1)  replacement  of  L  by  the  complex  L  =  L  —  2 (2)  a  finite  dispersion  of  u/n*; 
and,  if  a(z,  z‘)  is  known  in  detail  in  the  surface  region,  the  detailed  field  and  density 
distributions  in  the  surface  region. 

We  now  turn  to  resonant  response  to  a  uniform  external  field,  Eoe~,ut.  Eqs.  (7), 
(9)  and  (10)  remain  unchanged,  but  in  Eq.  (11)  E(z)  is  replaced  by  E(z)  —  E0  and 
Eq.  (12)  becomes 

jjf  \E(z)  -  £„]  =  £+'  dz'a(z,z'-,u)E(z').  (31) 

Near  u  =  wp  we  can  write  this  as 


^ E{z )  -  f  dz'^r^-EiA  +  r,  -±-E(z)  -  / dz,T-^^-E{A  = 

47T  J  t  47T  J  l  4x 


(32) 


where  (7(2,/)  =5  (7(2, 2';  wp),  7(2,2')  =  [8(7(2, 2';  u)/du]^r  and  7  =  a  -up.  On  the 
right-hand  side,  r\  can  be  neglected. 

Setting  Eq  =  0,  for  a  moment,  gives  us  the  equation  for  the  slightly  complex 
normal  modes,  En(z),  and  eigenvalues,  rjn  =  7^  +  it;".  The  kernels  <r(z,  z')/i,  T{z,z')/i 
are  symmetric  and  (slightly)  complex.  The  En(z)  satisfy  the  orthonormality  relations 
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L  +  c  rL+c 


cr. 


dzdz'E . 


47T 


£„(c')  =  6mn. 


(33) 


We  now  substitute  the  Ansatz 


n 


(34) 


into  (32),  giving 

■  £A„(,, -,„)£„(*)=££„.  (35) 

n 


Using  the  orthonormality  relation  Eq.  (33)  leads  to 

UpEfi  ri+c  r£+e 


An  = 


jLn  fL'+c  rL+c 

rb/«  */. 


4*(q  -  *?n) 


1  a  /x  r(*>*') 


(36) 


Only  even  modes,  n  =  2,4 . . .  are  excited.  A  well  defined  resonance  occurs  for  i\  as  r\'n 
if  t j"  <  }»7^  -  j/J li2).  In  this  case  as  single  term  in  (34)  is  dominant,  giving 


En{z)  =  Eq 


Up  [(?  -  Vn)  -  «>7n] 

4*  [fo  -  ‘ 


(37) 


In  order  to  resolve  the  discrete  modes  experimentally,  the  spacing  between  the  modes 
must  be  larger  than  the  width.  From  Eqs.  (22)  and  (23)  we  find  that  the  spacing 
and  width  are  given  by 


(»?«  ~  Vn-l) 


V 


it 

n 


(38) 

(39) 


Thus,  as  n  increases,  the  width  gains  on  the  spacing.  For  a  film  of  width  L'  =  30  A, 
d'fy  =  d"  =  1  A,  and  a '  =  0.57,  a"  —  0.025  (Al),  the  highest  well  defined  resonance  is 
n  as  20. 

We  plan  to  address  the  considerably  more  complex  theory  of  modes  with  finite 
wavevector  p  (Eq.  3)  in  the  near  future. 
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FIGURES 


FIG.  1.  Positive  background  density,  n+(r)  (dashed  line),  and  equilibrium  electron 
density,  n_(z)  (dotted  line). 


FIG.  2.  The  limiting  (w  — ►  up)  form,  E(z),  and  mode  density  n(z )  in  the  surface 
regions.  The  mode  density  n(z)  approaches  the  constant  value  r»o  in  the  bulk. 


FIG.  3.  (a)  The  real  part  of  the  electric  held  of  the  three  lowest  bulk-like  eigenmodes 

near  the  surface.  In  the  surface,  the  fields  all  approach  the  limiting  form  E(z).  Insert  The 
real  parts  of  the  three  lowest  bulk-like  eigenmodes  for  0  <  z  <  1/2.  In  the  bulk  region,  the 
fields  are  sinusoidal.  The  vertical  dashed  line  marks  the  position  of  a. 


FIG.  4.  (a)  Real  and  (b)  imaginary  parts  of  the  plasmon  eigenfrequencies.  The  straight 
lines  are  the  dispersions  Eqs.  (22)  and  (23)  with  V  —  34.8  A  and  dy,"  =  6  x  10~2  A  obtained 
from  Eq.  (20). 
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